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0. INTRODUCTION 


The Grant NGR-06-002-147, "Theoretical Investigations on Plasma 
Processes in the Kaufman Thruster," is concerned with I) the sputtering 
of the accelerating grid, II) the sputtering of the cathodes of the 
hollow cathode and neutralizer discharges, III) the deposition of the 
sputtered atoms on system components such as the solar energy collectors, 
and the IV) hollow cathode and neutralizer discharge characteristics. 

The progress made on these subjects in the period from 6.1.73 to 6.15.74 
is communicated herein. 

In Part I, an analysis of the sputtering of metal surfaces and grids 

3 

by ions of medium energies (e - 10 eV) is given. The sputtering is 
explained by discontinuous, nonlinear thermal waves (generated by the 
impinging ion) which produce a spatially concentrated emission of metal 
atoms under strong nonequilibrium conditions. It is shown that the con- 
ventional parabolic (approximation) heat conduction equation can not 
describe the transient transport of heat in metals at high temperatures 
(t £ 300°K) and has to be replaced by an exact, nonlinear, hyperbolic 
wave equation for the temperature field. This approach leads to a 
theoretical prediction of the threshold energy for sputtering and to a 
quantitative theory of the sputtering rate. As concrete applications, 
i) the number of atoms sputtered from the accelerating grid by a charge 
exchange ion beam and ii) the sputtering of system components by micro- 
meteorites are discussed briefly. 

In Part II, a quantum statistical and a perturbation theoretical 
analysis of surface sputtering by ions of low energy (e < 10 eV) is 
presented. Both approaches lead essentially to the same expression for 


1 


2 


sputtering rate, l.e. dependence on the ion energy, atom density of the 
solid, the atom and ion masses, and scattering cross section. The 
theoretical sputtering rate formula agrees well with experimental data, 
in particular as to the threshold energy and the energy dependence. As 
an application, the number of atoms sputtered from the cathode of low 
pressure discharges is calculated. The underlying model assumes that 
a quasi- thermal ion beam is formed in the potential drop of the cathode 
sheath. 

In Part III, the deposition of sputtered atoms on system components 
is treated. The transport model assumes that the sputtered atoms do 
not interact with themselves or any of the plasma particles (transport 
by free atomic flow) . Analytical formulae for the deposition rate 
are given in the case of uniform, nonuniform parabolic, and arbitrary 
nonuniform emission sources. Only such system surfaces are considered 
which can be seen along straight lines from the emitter. 

In Part IV, the theoretical efforts in determining the potential 
distribution and the particle velocity distributions in low pressure 
discharges, such as the hollow cathode and neutralizer discharges, 
are briefly discussed. Although two additional months were invested 
in the resolution of these problems, it was not possible to complete 
it because of mathematical difficulties. It is shown that the 
description of a collisionless electrical discharge leads to a 
nonlinear boundary-value problem for the coupled Vlasov equations 



and the Poisson equation for the electron and ion components, which 
has functional boundary conditions. In spite of a significant effort, 
it was not possible to determine the specific discontinuous functional 
solutions which satisfy the nonlinear functional boundary-value problem. 
The purpose of the investigation is to calculate the potential distri- 
bution, in particular the cathode and anode falls, and the electron and 
ion velocity distribution*.. The velocity distribution of the ions is 
of interest in connection, with the sputtering at the cathode. It is 
hoped that this investigation can be completed at a later date. 

The Investigations reported herein represent preliminary 
communications. An extended version of this work will be communicated 
in form of publications. In the past research period, the following 
investigations were published: 

1. H. E. Wilhelm, Transient Ion Neutralization by Electrons, 

J. Appl. Phys. 44, 4562 (1973). 

2. H. E. Wilhelm, Intercomponent Momentum Transport and 
Electrical Conductivity of Collisionless Plasma, Can. J. 

Phys. 51, 2468 (1973). 

3. H. E. Wilhelm, Nonlinear Theory of Electron Neutralization 
Waves in Ion Beams with Dissipation, Phys. Fluids 17 (1974). 
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X. CONTRIBUTIONS TO VOLUME SPU1TERING 
1. PROBLEM AND FOUNDATIONS 

In the evaluation of the sputtering of metal surfaces (cathodes, 

grids), two classical problems are encountered, i) the determination 

of the energy distribution of the sputtering ions at the metal surface 

(kinetic problem) and 11) the calculation of the number of atoms ejected 

by an ion of given energy based on a physical model for the sputtering 

mechanism. The phenomenological approach to the sputtering process by 
1 - 2 ) 

von Hlppel-Townes 7 assumes a Gaussian temperature distribution T(r,t) 

around the point of impact of the ion at the metal surface which reaches 

to infinity (infinite speed of heat propagation) and flattens out as 

time increases, T(r,t) 0 for t ». The vapor pressure P(r,t) of 

1 - 2 ) 

the metal is assumed to adjust Itself instantaneously to this 

transient temperature distribution in accordance with statistical 
equilibrium mechanics, 

P(T) = (18rMu) 2 ) 3/2 /kT)“ 1/2 exp(-E g /kT) 

where E g is the sublimation energy and w D = k0 D /2h a frequency 

3) 

related to the Debye-temperature 0 — 7 . This approach is unrealistic 

in the treatment or the thermal dissipation of the ion energy and 

assumes a physically unrealizable transient metal vapor equilibrium. 

As one sees from the above formula, the von Hippel-Townes model does 

not give a threshold energy for sputtering (since P > 0 for any T > 0) 

4-5) 

as observed in experiments. 

The other theoretical approaches are based on considerations of 

6 o') 

momentum conservation ("focusing, collision sequences"). 7 They 
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predict a threshold energy for spattering which depends strongly on the 
masses of the ion and the metal atom.— — ' This result is in direct 
disagreement with the experimentally determined threshold energies 
which are of the order of the dislocation energy of an atom in the 
metal lattice.—”—^ 

1—2 6 — 8 ) 

The theoretical models used in these attempts^ — 1 at 

explaining the sputtering process contain phenomenological parameters. 

Their success in explaining experimental observations appears to be 

due mainly to a proper adjustment of the phenomenological constants 

in each case. In the following, we try to develop a volume sputtering 
2 

theory (e > 10 eV) which is free from phenomenological parameters. We 

show that the impinging ion, which penetrates through a certain number 

of atomic layers in the metal, generates a non-linear, discontinuous 

thermal wave. As a result of the high concentration of energy behind 

the wave front, the thermal wave produces a mechanism which breaks the 

atoms out of their bound places in the lattice. The metal atoms in the 

volume overrun by the thermal wave are emitted until the energy behind 

the wave front has decreased down, to the dislocation energy of an atom. 

This concept leads directly to the correct threshold energy for 

2 

sputtering. A theory of surface sputtering (e < 10 eV) is presented 
in Section II. 

In the slowing down process of an ion penetrating into a metal, 

its kinetic energy e is dissipated nearly homogeneously along its 
9) 

path of length L.— Accordingly, the energy expended per unit path 
length is in this approximation 


6 


e ■ e/L [erg cm 

where c is Che known ion energy, while L can be calculated from 

9 ) 

slowing down the«J>*y for charged particles in solids.— The energy e 

of the ion appears quasi-lnstantaneously in form of thermal energy due 

-19 

to the high density of the metal (relaxation time - 10 sec). Titus, 
a cylindrical thermal wave is generated by the ion with the path L 
as symmetry axis. (In cuse of low ion energies, when only one or a few 
atomic layers are penetrated by the ion, an essentially semi-spherical 
thermal wave is generated around the point of Impact.) 

The transport of heat in a metal is described by the relaxation 
equation for the heat flux q and the conservation equation for the 
thermal energy density pcT (p « mass density, c = specific heat, 

T = absolute temperature) . These equations are derived as moments of 
the Boltzmann equation, and are: 

f£ - - 7 q - ~ ro > (1) 

pc “ = -v.q , (2) 

where 

x — t(T) 5 x o T m , [r Q ] = sec deg - ™ , (3) 

X “ X(T) = X Q T n , [X Q ] = erg cm ’'’sec ^deg n ^ , (4) 

are the relaxation time of the heat flux and the thermal conductivity of 

the metal, respectively, The temperature dependence of r and X can be 

modelled in wide temperature ranges by simple power relations (m, n * 0) . 

+1 

Theory and measurements indicate that X - T at low temperatures, 
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*»2 Q 

A - T at Intermediate temperatures* and X - T at high 

temperatures.—^ Since the electrical conductivity of the metal, 
2 

a " ( n e e / m e ) T *» is proportional to the momentum relaxation time, 
T io given in termB of X by the Wiedemann-Franz relation—^ 


t(T) 



( 5 ) 


where the electron density n e is to be considered a constant 

(pc ■ p q c o for a quasi-incompressible metal). Consideration will be 

given exclusively to cylindrical thermal waves. 



8 


2. PARACOLIC THERMAL WAVES 

In contemporary heat transport theory,^ It is standard to 

assume that a temperature gradient VT produces instantaneously a heat 

flux q, l.e. Eq. (1) is replaced by q ■ ~XVT. Combining this relation 

11 ) 

with Eq. (2) gives the usual parabolic heat conduction equation. — 

The initial-value problem for cylindrical thermal waves becomes in 
this approximation: 


II 

3t 


„ I J_ / T n v 8T. 

a r 3r (1 8r' 


( 6 ) 


where 


2 nj T(r,t) rdr ■ Q 
0 


CD 


expresses the conservation of energy deposited per unit length by 
the ion» and 


a 2 X Q /pc 
Q = e/pc 


[cm 2 sec’*'*' deg n ] 


[deg cm ] 


( 8 ) 

(9) 


Eq. (7) is mathematically equivalent to the initial condition, 

T(r,t « 0) = Q 6(r)/2irr. The parameters a,Q and x,t permit the 
formation of a single nondlmensional combination, 

i 

i 

£ - r/(aQ n t) 2(ltHL) . (10) 


which has the meaning of the similarity variable of Eqs. (6) - (7). 

For dimensional reasons, one makes for the temperature field the ansatx 


1 

T - (Q/at) n+1 f(|) 


(ID 
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Eq. (11) reduces Eqs. (6) - (7) to an ordinary nonlinear problem for 
the nondlmenslonal function £(£)» 


where 


2 0+1) A <f n t If) + C 2 § + 2Et - 0 


2 it/ f(0 5 d? v 1 
0 


(12) 


(13) 


Eq. (13) has a closed form solution which Is dlscpntlnuous [H(x) *» 1, 
x > + 0; H(x) ■ 0, x <. - 0] : 

1 

■ (jTSoT «« - «V 11 «<T 5 > < w > 

where 


-2<n+l) 


i i 


l dn 


l.e. 


i n 

2 ,.n, .n+1 n+1 

- (4 M — 


(15) 


by Eq. (13). It Is physically more illustrative to rewrite the 
temperature distribution of the thermal wave as 


A 2 — 

T(r,t) - T(t)[l - — • — ] n H[R(t) - rl 
R 2 (t) 


(16) 


where 


T(t) - 7 r[n/ 4 (n+l)] n 5 2(n+l)/n T/t) 


T(t.) = Q/nR^t) 


, (17) 

, (18) 


and 
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R(t) - 5 o (aq n t) 2 ^ n+1) . (19) 

Eq. (16) indicates that the temperature in the wave drops 
discontinuously to zero at r ® R(t), the position of the wave front. 
A nearly homogeneous concentration of thermal energy exists behind 
urte thermal wave front, 0 <_ r £ R(t), which advances with the speed 


dR(t) 

dt 


2 (n+1) 


t _1 (aq n t) 


1 

2 (n+1) 


( 20 ) 


It should be noted that the parabolic solution in Eq. (16) 
diverges for n < 0 at r ** R(t) , and reduces to an unrealistic 
Gaussian distribution which extends to infinity at any time 0 < t < » 
in the linear case, n = 0. Accordingly, for media with \ decreasing 
with increasing T or constant X(n £0), Eq. (16) does not represent 
a useful approximation. As will be shown, these difficulties are 
removable by means of Eq. (1) which takes into consideration the 
physically required relaxation in any transient heat flow. 

The propagation of a parabolic thermal wave in barium oxyde (BaO) 


is shown in Fig . 1 (n = 3 ) . 
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3. HYPERBOLIC THERMAL WAVES 

In n transient heat transport process, a temperature gradient 
produces a heat flow after a finite relaxation period. In accordance 
with Eqs. (1) - (3), the cylindrical thermal wave generated by an ion 
of energy e around its path L in a metal is described by the 
hyperbolic initial-value problem; 


3£ n _ 1 _ X. 3T 

Dt T ^ t Jr 


( 21 ) 


3T 1 3 , . 

pC 9F “ “ 7 97 (rcl) 


( 22 ) 


where 

m 

2iTpc / T(r,t) rdr = e . (23) 

0 

Sputtering is produced in metals exclusively at temperatures T > 273°K. 
In thlB so-called high temperature region, X is constant whereas x 
is inversely proportional to T [Eq. (5)], i.e. 10 

x - t o T _1 , [t q ] = deg se , (24) 

X = X Q T° , [X o ] = erg cm 3 sec 3 deg ^ . (25) 

Since [e/X Q T o ] = 0 and [X Q /pc] = cm sec , a (nondimenslonal) 
similarity variable results from the dimensional parameters e", X q , 
t o> pc and r,t in the form 

E = r/[ (X Q /pc)t] :l '^ . (26) 

For dimensional reasons ( [t q 3 = deg sec, ((X pc) 1 ^ 2 r ] = erg cm” 3 cm 
-1 3/2 

sec sec ), the temperature and heat flux fields are subject to 
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the self similar transformations, 

T(r,t) » T o t _1 f(5) , (27) 

q(r,t) - U Q pc) 1/2 r o t” 3/? * g«) . (28) 

These equations reduce Eqs. (21) - (23) to a problem for ordinary 
nonlinear differential equations: 

5 ff + 3g = 2f[g + ||] , (29) 

Id - l_df 

£ d5 " f + 2 S d? , (30) 

where 

2tt / f($) K d£ = e, e = e/1 t , (31) 

0 

and f(£) and g(£) are nondlmensional, Eq. (30) is readily 
integrated , 

a Z t - 2£g) - 0, i.e.: 5 2 f-2£g = C Q , 

whence 

g - | Cf , (32) 

since C Q =» 0 by the condition g = 0 for £ * 0[q(r ss 0,t) = 0 for 
reasons of symmetry]. Elimination of g from Eq. (29) by Eq. (32) 
yields 

~ 4f) Jf = 2?f(f-2) , (33) 

or 

a 2 - 4f) = f (f - 2) (34) 

dr 
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By multiplication with the integrating factor, 

M (f ) « f“ 1/2 (f-2)" 3/Z , (35) 

Eq. (34) is transformed into the complete differential, 

2 1/2 

dS(g 2 ,f) ° . df + — ' 1/2 dg 2 = 0, (36) 

f A// (f-2) A/Z (f-2) A/ ^ 

Accordingly, the solution f(£) is given implicitly, by the Integral 
S(5 2 ,f) - C of Eq. (36): 

5 2 - 8 + [(f-2)/f] 1/2 {C - 8£n[f 1/2 +(f-2) 1/2 ]} . (37) 

This result yields directly 5 “ 5(f) and by inversion the solution 
f 55 f(5), which is symmetrical, f(+5) = f(-£). In particular, Eq. (37) 
indicates that a real solution f(5) > 0 exists only in the interval 

for 8>? 2 >0 , (38) 

where 

f = f min - 2 for 5 - + 2 /2 , (39) 

f = f = f for 5=0 , (40) 

whereas 

f = 0 for |5| > 2 /Z , (41) 

by Eq. (34), l.e. f(5) is discontinuous at % - + 2 ^2 . The value 
f Q = f(5=0) is related through Eq. (37) to the integration constant C, 

C/8 = + (f Q - 2) 1/2 ] - [f o /(f Q - 2)] 1/2 . (42) 

The energy conservation relation in Eq. (31), which determines C and 


w ... 
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thus f Q , becomes 
2 2 

/ f ^jL df = e/ir , (43) 

f o 

2 

by change of the Integration variable. Substitution of d£J /df and 
C in accordance with Eqs. (37) and (42) leads to the transcendental 
equation, 

•Mf^ 2 + (f Q - 2) 1/2 ] Un/2 + [f Q /(f o - 2) 1/2 ] - | £n[fj /2 + (f Q - 2) 1/2 ]} 

+ -(f o - 2) - { [f 0 /(f Q “ 2)] 1/2 + | ln/2 } ln/2 - e/l6Tr , (44) 

which gives f as the first real root f >2. The left side L(f ) 
o o o 

of Eq. (44) assumes the value L = 0 for f ** 2 and Increases 

monotonically with increasing f > 2 so that L » for f 

0 o 

Hence, a real root 2. < f < " exists for any given 0 < e < » . 

It should be noted that the function f(£[) loses its uniqueness 
at sufficiently large e ^ e . According to Eq. (37) 

2 2/2 

[df/dr], ^ = — , 0 < e « 1 . (45) 

f = 2 + e C-8 ln/2 

\ 

Hence 

[df/d£ 2 ] £ =2+e <0 for C < 8 ln/2 . (46) 

n 

This means that f(£) is i) a unique function of $ for C < 8 Zn/2 

but ii) a multivalued function of £ 2 for C > 8 tn/2. Since C = C(f ) 

o 

[by Eq. (42 ) 3 and f Q = f Q (e) = £ Q (e/\ Q T Q ) tby Eq. (44)] increase with 
increasing f Q and e ^ e, respectively, multivalued thermal flow 


• T 

J 
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appears for e > £ where 

U[f 0 (&)] - 8 InSl . (47) 

A 

Accordingly, for an energy release e > e, where the critical value 

A 

e = & X o t q is determined by Eqs, (42), (44), and (47), the energy e* 
is no longer propagated through an ordinary nonlinear thermal wave but 

through a thermal shock wave. As expected, thermal shock waves 

___ * 
occur at energies e above a critical value e. 

By combining the analytical solution f = f(£) in Eq. (37) with 

Eqs. (27) - (28) and (32), one obtains the fields T(r,t) - t 2 f(£) 

and q(r,t) ~ t~ 3 ^ 2 £f(£) where ? - rt -1 ^ 2 by Eq. (26). T(r,t) and 

q(r,t) are decreasing with increasing t at any fixed point 

0 < r < r(t) within the thermal wave. These fields are discontinuous 

at the wave front, 

R(t) = 2/2 U o /pc) 1/2 t 1/2 

by Eqs. (20) and (39). The speed of the wave front is 
= a a /p C ) 1/2 t' 1/2 

at o 

Accordingly, the wave spreads out radially with time at a speed 
decreasing with time. It is interesting that R(t) and dR(t)/dt are 
independent of e and t q for the particular T - dependence of X 
and t in Eqs. (24) and (25). In this case, s and t affect only 
the height of the distribution f(£), i.e. f increases with increasing 
e and decreasing t q [Eq. (44)]. 


(48) 


(49) 
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The propagation of a hyperbolic thermal wave in wolfram (W) Is 
shown in Fig. 2 (m 55 -1, n ■ 0) . It is seen that the hyperbolic wave 
exhibits an extremely steep wave front in comparison to the parabolic 
wave (Fig. 1). 


= 2.4 x 10 
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4. APPLICATION TO VOLUME SPUTTERING BY IONS 

For materials which have a thermal conductivity \ T n , n > 0, 
such as metal oxides (e.g. BaO) , glasses, graphite, etc., at high 
temperatures, the parabolic equation gives an approximate description 
of thermal waves. The parabolic solution [(Eq, 116)] does not indicate 
the occurrence of shock waves at high ion energies e and its approxi- 
mate validity is, therefore, questionable in this energy region. In 
the case of pure metals, which have a thermal conductivity X 'v T^ and 
a thermal relaxation time t * T~\ only the hyperbolic system in 
Eqs. (21) - (22) provides a physically acceptable description of 
thermal waves. The hyperbolic solution [Eq. (37)] is applicable for 

_ _A_ 

ion energies e < c [Eq. (47)], since it becomes multivalued for 

— A 

e > e (shock waves). 

A thermal wave of cylindrical symmetry represents a first 
approximation to the actual thermal waves produced by sputtering ions 
in materials. Deviations from the cylinder symmetry are due to end- 
effects at the point of impact and the end of the ion path (in particu- 
lar at low ion energies e), nonuniform slowing down, anisotropies in 
the directions of the most dense atom arrangement. In these directions, 
the probability for momentum transfer is largest so that the resulting 
crater in the material resembles more a cone than a cylinder. 

The number of atoms emitted by an ion of given energy e is 
mainly determined by the energy conservation equation (whereas the 
calculation of the spatial distribution of the expelled atoms would 
require in addition consideration of many-body momentum conservation 



20 


In the solid). For this reason, the number of atoms sputtered by an 
Ion may be evaluated macroscoplcally by means of the thermal wave 
concept. In this approach, the collective many-body Interactions which 
produce microscopically the thermal energy transfer are contained in 
the thermal conductivity \ and the relaxation time t. The crystal 
bonds of the atoms behind the wave front R are broken so that the 
atoms are emitted as long as the average particle energy McT(t) in 
the thermal wave is larger than the effective threshold energy E q 
for sputtering, 

R(t) , ~ - 

McT(t) = 2irMc / T(r,t) r dr/flR (t) >_ E q , E q « aE Q . (50) 

where M = p/N is the mass of an atom in the solid. The correction 

factor a takes into consideration that on the average a fraction of 

the energy e goes into kinetic energy of the expelled atoms, i.e. in 

general 1 < a < 2. According to experiments, the true sputtering 

threshold E is a material constant which is independent of the mass 
o 

4 - 5 } 

ratio of the atom and ion.- The sublimation energy E. is the 

b 

energy required on the average for the removal of an atom from the 

surface of a polycrystalline solid, E g - < E^ 1 * > ijk (ijk designates 

the orientation of the surface) . If the atom is expelled from within 

the solid, then E Q is equal to the dislocation energy E d of an 

atom. E„ is the energy required for i) the removal of an atom from 
d 

its position in the lattice (~2E g ) and it) its stable transfer to an 
interstitial lattice position (~2E g ) , i.e. the threshold energy E q 
for sputtering is proportional to E g , 
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E - HE , h - 4 . (51) 

o s 

In Table I, the experimentally observed threshold energies E q for 

4 - 5 ) 

some technically interesting metals— — are compared with the theoreti- 
cal values following from Eq. (51) for h " 4 [note that E q (experi- 
mental) was obtained by independent logarithmic extrapolation of the 
data in references 4-5]. The discrepancies between the experimental 
and theoretical values of E q lie well within the experimental 
uncertainties. 

TABLE I. Experimental and Theoretical Threshold Energies. 


Metal 

E q ( experimental ) 
eV 

E s 

eV 

E q (theoretical) , h*4 
eV 

A1 

16 

3.3 

13.2 

Cu 

20 

3.5 

14.0 

Mo 

26 

6.2 

24.8 

T.i 

35 

8.0 

32.0 

W 

36 

8.8 

35.2 


Equation (50) defines a maximum wave front radius within which the 
crystal bonds of the solid are broken. These atoms are expelled with a 
relaxation time of the order t £ a tf/E Q . At this phase of the expansion, 
the thermal wave collapses as its energy e has been consumed in 
expelling the atoms. According to Eq. (50), R(t) is given by 

R(t) = [e/irE o (p/M) ]"* (52) 

for both the parabolic and hyperbolic thermal waves. The number Z of 
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atoms sputtered by an ion of energy c is on the average (a) 

Z - nR 2 (t)L(p/M)H(e - E 0 ) - <e/£ 0 )H(e - E Q ) , (53) 

H(e - E q > - 1, e > E q 

* 0. e < E . 

o 


Equation (53) is based on a discontinuous propagation of the ion energy 
e which is provided by the discontinuous thermal wave in the continuum 
picture (e » E q , Z » 1). Substitution of Eq, (52) into Eqs, (19) 
and (48) yields the time it takes the parabolic (p) and hyperbolic (h) 
waves to propagate to the critical radius R(t), 


~r 1 / n y n+1 lie. n+1 e 

S ■ 4? W <-p r o 


B / X n E n 
O O 


n+1 


t. ■* T" “ “ g/X E 

h Sit - X o o 

E o 
o 


(54) 

(55) 


In applications to ion sputtering, it is to be noted that Eq. (53) 
is valid for not too large ion energies. At high ion energies, the 
ion penetrates so far into the solid that only relatively few volume 
atoms are emitted. 

Sputtering by ions with energies significantly larger than the 
threshold energy (e » E q ) occurs at the accelerating grid of ion 
propulsion devices.— —^ The velocity distribution of these ions 
at the surface of the accelerating grid (z = 0, 0 < r £ R) may be 
simulated by a Gaussian of the form 

t in ^ ..3/2 -m(v - <v>) /2<e> te . fS 

f (v) = n(m/2ir<e>) e (56) 


r 


! 
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where 

n - n(r) ■ Ion density at z“0, O^r^R 

<c> " <e(r)> ■ average random energy at z«0, O^r^R 

<v> ■ <v(r)> « directed beam velocity at z«»0, O^r^R 

*y 

v ■ individual ion velocities relative to grid 

If che beam is in thermal equilibrium in its center of mass system then 

one has <e> - lcT(r)« The number of grid atoms sputtered per incident 

2 

ion of energy e “ h mv is by Eq. (53). 

Z(h mv^) « [h mv^/E Q ] H(*i mv^ - E q ) . (57) 

Let 0 designate the angle between the vectors v and <v>. The 
number of atoms expelled per unit area of the metal surface and unit 
time is at the radial location 0 < r < R 


dN _ / m k3/2 m 


2E 


2n -nil 00 „ ,-y *.2. k 

rtf 2 -m(v~<v>) /2 < e>, 2 

[ ) J v e (v cos6)v sln9 d+dOdv 

° , (58) 

h 


0 0 v 


v 0 5 (2B 0 /«) 


After the trivial ^-integration and the substitutions, 


cos8 = T, 
/g v = x, 


dt = -sin0d0 
dx » dv 


( 59 ) 


24 


and 


3 - m/2<e>, 
x o - /B <2E /m)* 2 , 


y = |<v> 


y 0 " x 0 “ Y 


(60) 


Eq. (58) becomes 


|| - n(m/2ir<e>) 3/2 <mn/i o )0" 3 J 


(61) 


where 


1 » , 2.2 
j . / / e - <x + v 

0 x_ 


2YXT) TdT X 5 dx 


(62) 


whence 


2 oo 2 

_ , -2 -y ,* 3 -x 

J - -$y e J x e dx 


+ % 1 / (y + Y) 4 e y dy 

y o 


- J 4Y~ 2 7 (y + Y) 3 e~ y 


dy 


(63) 


i.e. 


. 1 -2 ... , 2. -(x 2 + y 2 ) 

J = g Y (1 + X' 0 ) e o 

+ ^ Y(5 + 2y 2 ) [1 - *(y 0 )] 


-|y" 2 [(1 " 5y 2 - Sy 4 ) -12y 3 y o 


2 2 3 , -v 

+ (1 - 8y )y„ “ 2-yy ] e y o 


( 64 ) 
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Since 

Y » (.h m<v> 2 /<e>)' i , x Q “ (E Q /<e>)'* , 

y o” <V <e> ^ - Qi m<v> 2 /<e>)' 5 , (65) 


one obtains for the number of atoms sputtered per unit area and unit 
time [Eq. (61)] explicitly: 


ft - »<* ^> h 


(66) 


where 


E . -<^ 2 3/2 


E 


a + £> °' (E ° + k ” <v> )/<e> 


+ m<v> /<£>) (5 + 2Qi m<v> ■) 2 ] } 


<f>> 


+ 2 i 2 2 2 3/2 

- {1 - 5 0s m<v> /<e>) -8 (h m<v> /<e>) - 12 (h m<v> /<s>) 

♦ [ (E_/<e>)^ - (h m<v> / <e>)' 2 3 1 + [1 - 8(% m<v> /<e>) ]((E /<e>)' 5 ~ 


(•| utfv? 2 /^)**] 2 - 2(^5 m<v> 2 /<e>)' 5 [(E o /<e>)^ 


y 1 2 * v ^5-, 3 >. 

0i m<v> /<e>) ] } 


-[ (E /<e>) 
e 0 


^ 2 

(% m<v> /<e>) 



(67) 


This result indicates that dN/dt varies in a rather complicated 

** 

way with increasing <e>, and <v> for fixed E q . In Fig. 3, 

H(E Q ,<v>,<e>) is shown quantitatively in dependence of the energy 
2 2 

ratios y and x q [Eq. (65)]. With exception of the region of low 

beam velocities <v> [compared to the average random velocity 
U 2 -1 

(2 <e>/m)^] , i.e. y < 10 , H(y.x q ) increases with increasing 

2 + 2 2 

Y - <v> /e and decreasing x q - E Q /e. 
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5. APPLICATION TO VOLUME SPUTTERING RY MICROMETEORS 


The nonlinear thermal wave theory presented explains also the 
sputtering produced by micrometoors impacting on system components of 
ion rockets and space ships in general. The above considerations on 
ion sputtering are transferable to sputtering by micrometeors . In the 
case of an Impinging micrometeor, the energy e * eL [Eqs* (7) and 
(23)] has to be interpreted ns 

e = h Miiv** , (68) 


where v is the velocity of the center of mass of the micrometeor, W 
the number of atoms it is made up of and m the average mass, 


m 


1 

N 



W m 
s s 


(69) 


W is the number of atoms of mass m contained in the micrometeor. 
s s 

The main chemical components of micrometeors are iron (Fe) or stone 

(Si0 2 ) ^ Their speeds v range from 10 km sec” 1 to 10^ km sec” 1 .— ^ 

The radii of micrometeors range from r . = 10~^ cm to r 5 10” 1 cm, 

min max * 

and their mass per volume of space greatly exceeds that of all other 
18 ) "1 

meteors — ' (r > 10 cm). For this reason, micrometeors are most likely 
to hit ion rockets and space ships outside of the atmosphere of the 
planets. For micrometeors, the effective sputtering threshold is 
about [Eqs. (50) - (51)] 

E = aE = ct4E s 8E , (70) 

o o s s * V ' 

i.e. eight times the sublimation energy E . The sublimation energy is 

s 

E /in = 7 x 10 10 erg/gr for Fe and E /m = 14 x lO 1 ^ erg/gr (including 
s s 
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dissociation energy) for SiO^* According to Eq. (70), the critical 
minimum speed for n micrometeor to vaporize itself on impact is 

v min " <2 *■'»* ■ <71) 

Hence, v. a (2 x 8 x 7 x = 10^ cm sec -1 = 10 km sec ^ for a 

min 

pure Fe-mlcrometeor , and v (nin a (2 x 8 x 14 x 10*°) ' 5 = 1.5 x 10 6 cm sec 
** 15 km sec ^ for a pure SiO^-micrometeor. In order that the thermal 
wave vaporizes also a significant volume of the material of the system 
on which the micrometeor impacted, the speed of the latter must 
satisfy the basic inequality 


v » v 


min 


(72) 


As the above examples indicate, this condition is satisfied for a 

2 - 

significant percentage of micrometeors which have speeds v = 10 ton sec 

and larger (note that the sublimation energies E of all solid 

s 

materials are of the same magnitude-of -order, i.e. a few electron volts). 
The number of atoms Z sputtered by a micrometeor with a speed 
v >> v m ^ n out of the target material is [Eq. (53)] 

Z = e/E Q - JiWmv 2 /E o . (73) 

-4 8 -1 
For a Fe-micrometeor of radius r - 10 cm and speed v = 10 cm sec , 

one has e /3)mv 2 /A = 2 x 10 ^ 9.3 x 10 10^/10 ^ s 2 x 10^ 

-23 3 

erg (atomic volume A £ 10 cm ) . On a wolfram target (E o = 35 eV, 

-12 -10 

E = 8 x 35 x 1.6 x. 10 5 3 x 10 erg), this micrometeor would 

sputter Z = e/E Q ® 2 x 10"*/3 x 10 ^ 5 10^ atoms. The corresponding 

— 1/3 

sputtering crater has an extension of the order R = (ZA) = 

nn 15 . -23.1/3 . ln -3 

(10 x 10 ) = 2 x 10 cm. 
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The formulas derived for thermal waves and volume sputtering 
in sections I, 3-5 are similarly applicable to the unloading and 
sputtering by micrometeors. The necessary modifications are defined 
in Eqs. (68) - (73). 
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II. CONTRIBUTIONS TO SURFACE SPUTTERING 

At sufficiently low energies of the incident ions, exclusively 

surface atoms of the solid are sputtered. The experimental data on 

sputtering of metal surfaces indicate that the average number S(E) of 

15 

atoms sputtered per Incident ion of energy E lies in the interval-” 

0 <_ S(E) <1, E q < E < 10 2 eV , 

where E q is the threshold energy for sputtering of surface atoms. 

The measured S(E) - curves can be fitted by analytical expressions 
25 

of the form— 7 

S(E) = a(E - E Q ) n , a = const, n = 2 , 

2 

at low energies, E q <_ E < 10 eV. It is shown that this simple 
sputtering formula can be explained theoretically by means of a 3-body 
sputtering mechanism involving the ion and two surface atoms of the 
solid. By means of a statistical analysis and a quantum mechanical 
perturbation theory one finds independently that n=2 and that "a" 
is a weak function of energy E which can be taken to be a constant 
for E * E q , i.e. a(E) = a(E Q ). 

An ordinary binary collision between a surface atom of the solid 
and an ion incident normal to the surface can evidently not lead to 
sputtering since the atom does not acquire a momentum component In 
the direction of the external normal of the surface. Similarly, 
sputtering is not likely to occur for smaller angles of ion incidence if 
its energy is not large compared to the threshold energy for sputtering. 
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It Is evident that sputtering, at ion energies of the order of the 

threshold energy, is a 3-body process involving one ion and two 

surface atoms of the solid. At higher ion energies, however, sputtering 

will result mainly from higher order many-body interactions. 

By restricting the theoretical considerations to ion energies E 

2 

of the order of the threshold energy E Q , E Q 5 . E < 10 eV, sputtering 
is regarded as the result of an ion-atom-atom interaction. Furthermore, 
it is assumed that the solid is polycrystalline and has a sublimination 
energy which is on the average E fl - < E g iiik > where the average is 
taken over the randomly distributed surfaces (ijk) of the crystallites. 
In this case, the sublimation energy E g represents the average binding 
energy of a surface atom. In the 3-body sputtering process, the ion 
transfers i) the energy E to the atom which is expelled and 
li) the energy 2E or 4E to the other atom depending on whether 
the latter is pushed to an unstable or stable interstitial lattice 
position, as well as iii) kinetic energy. Accordingly, the average 
sputtering threshold should be 

E = J*(E +2E;+E + 4E ) = 4E o . (1) 

o s s s s' s 

In experiments which cannot detect individual but only a large number 

of sputtering events (e.g. sputtering of glow discharge cathodes and 

accelerating grids) , the threshold E q represents always an average 

value, i.e. not the absolute smallest possible binding energy of a 

surface atom which can be as small as E J^ k /5 where i'j'k* 

s 

designates that surface which has the smallest sublimation energy. It 
is interesting that the (average) threshold for surface sputtprlng is 
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equal to the threshold for volume sputtering since the energy for 
stable displacement of an atom within the solid is also 4E . 

When an ion of low energy as defined above hltB the surface of 
a solid, one of the following processes may occur: 1) the ion is 

reflected without energy loss by the bound surface atom it encounters} 
2) the ion collides with a surface atom and quasi-simultaneously with 
a second atom so that 3-body sputtering results. (More precisely, 
the designation "atom” should be used for the incident "ion" since the 
latter certainly recombines with an electron as it approaches the sur- 
face of the solid.) The total probability for the ion to Interact in 
either of the two ways with tie solid is 

P N - N 2/3 a(E) (2) 


where N is the number density of atoms in the solid and cr(E) is 
the (energy dependent) cross section for ion-atom scattering, let 
W^(E) and W 2 (E) be the probabilities for the processes 1) and 2), 
respectively. The relative probability with which sputtering occurs 
is then 


"s<» 


W^(E) + W 2 (E) 


W,(E) 

W^CE> > V E > « 


(3) 


Combining of Eqs, (2) and (3) yields for the sputtering rate, i.e. 
the number of atoms expelled on the average by one ion of energy E 
from the solid, 

S (E) = a(E)N 2 ^ 3 W (E) . (4) 

For the evaluation of W (E) , two methods, which are based on 

s 

different approximations, will be used. The total cross section 


I 
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o(E) is assumed to be known either theoretically or experimentally.— —^ 
One of the methods holds in the cose that the interaction by the ion 
can be treated qb a perturbation, whereas the other method holds for 
arbitrary strong interactions but assumes quantum-statistical equilib- 
rium among the final sputtering states. The latter assumption appears 
to be questionable at first sight since the small energy region under 
consideration permits only a relative small number of final states. 

As a justification it is noted that these different approaches lead 
essentially to the same result. 
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1. STATISTICAL ANALYSIS 

The basic assumption is made that the final sputtering states 
correspond to a quantum statistical equilibrium brought about by the 
strong Interaction between the ion and the surface atoms « In the 
processes 1) or 2), the ion Interacts with the surface of the solid 
within an area of the extension of the de Broglie wavelength, 

1 =» fi//2mE. For this reason, the spatial part of the phase space is 
taken to be-^ 

V » —■ R 3 , R » If/^mE) 15 . (5) 


In quantum-statistical equilibrium, the probability for transition 
into a final state 1 b proportional to i) the probability that the 
Interacting particles are simultaneously within V and ll) the 
density of final states dp/dE per unit energy. For a state con- 

-K *4* 

tainlng n independent particles with momenta p^, p 2 , •••> P n » w 


and dp/dE are given by 


/V.n dp _ r g nn d<i>(E) 
w ~ W > dE ~ i J hr 




dE 


( 6 ) 


n designates the normalization volume, fi > V, and 4>(E) is the 
volume of momentum space corresponding to the total energy E. 
Accordingly, the probability for transition into the final state n 
under consideration is 

W(E) = [V/(2TTfi) 3 ] n . (7) 


This equation represents the basis for the determination of the pro- 
cess probabilities W^E) and W 2 (E). Because of the conservation 
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1 awe, the n particles In the Interactions 1) and 2) are not 
Independent. This requires certain modifications of Eq. (7) which 
are explained in the applications below. 

W.(E) is defined ns the probability for the ion to be reflected 
at the surface of the solid without energy loss. In the center of mass 
system (solid) , the ion momentum is p * /2mE in the final state and 
the momentum space volume is <t>(E) - 4rp /3. According to Eq. (7), 
the probability for reflection is (n*l) 

W X (E) =■ [V/^irti) 3 ] 1 4ir/2m 3/2 E 1/2 . (8) 


W 2 (E) is defined as the probability for the 3-body sputtering 
state. In the center of mass system, the momenta of the ion (1), 
the sputtered atom (s) , and the second atom (a) can be chosen as 

, (9) 


-> -y -y 1 *»■ -v 
p 4 58 P> P Q T P “ 


1+.+ 
- 2 p + q 


so that momentum is conserved = Since the energy E q [Eq. (1)] 

is expended in the sputtering interaction, the total kinetic energy 
of the three particles is 


„ „ , 1.1 s -»-2 , 1+2 

E* = E - E q * (^ + p)p + „ 9 


M 


( 10 ) 


Equation (10) represents an ellipsoid with the axes sections 
{ [mM/ (m + 2M)]E*}^ 2 and (ME*) 1//2 in the six-dimensional space of the 
vectors p and q. Hence, the volume of the momentum space is 


0(E) 


■ —) 3/2 E * 3 

6 '“m + 2M ; L 


(ID 


37 


Owing to conservation of momentum, only two of the particle momenta arc 
independent, i.e, n“2. Substitution of Eq. (11) into Eq. (7) gives 
for the probability of the sputtering state 

W 2 a) - [V/ <2irR) 3 ] 2 4» 3 (^p^) 3/2 E* 2 . (12) 


With the assumption W 1 (E) » W 2 (E), one obtains from Eqs. (5) 
(8), and (12) for the relative sputtering probability the approximate 
expression 

2 

„ ; X. f - (M/nQ 2 _ ,3/2 _ (E _7 -V- / 13 n 

W s (E) “ 27 [ 1 + 2(M/m) ] J. (13) 


where V has been eliminated in accordance with Eq. (5). 
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2. PERTURBATION THEORY 


If the interaction energy between the ion and the surface atoms 
is sufficiently small, 1) the reflection of the ion by a surface atom 
without energy loss, and 2) the 3-body ion-atom-atom sputtering with 
an energy loss E > may be treated by perturbation theory. The 
probability rates for these processes are most conveniently determined 
by means of Fermi's Golden Rule.~^ 

In the case of the ion reflection, the magnitude of the momentum 
after the collision is p ** /2mE in the center of mass system. The 
probability per unit time of the reflection transition Is 




2tt | m (1) 1 2 do. 
1 "if 1 dE 


(14) 


where 

= /// ¥ f * H d 3 r 

H 


and 

% = [ft/(2irtl) 3 ] 2Tt*^m 3/2 E l/2 (16) 

dE 

are the matrix element of the perturbation (operator) in the 

Hamiltonian of the ion-atom system which causes the transition i->f 
and the density of final states per unit energy, respectively. 

and i|» f are the wave functions of the total system before and 
after the transition which are normalized for th-j volume £2. 

According to Eqs. (14) - (16), the probability for the ion reflections 
is per unit time 

Wl (E) = ^ [£2/(2TrK) 3 ] 1 |M^ ) | 2 4ir/2m 3/2 E 1/2 . 


(17) 
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The final state of the sputtering interaction consists of an ion 

1 1 ) r ^ 

of momentum p^ « p, a sputtered atom of momentum p s = - ‘jj p “ 9> 
and a displaced surface atom of momentum p a * - ^ P + 9 the center 

of mass system. Since 2^ P^ » ft, only two particle momenta are 
Independent which determine the statistics. The probability per unit 
time of the sputtering transition is 

V E > -f l M i £ 2> | 2 dE < 18 > 


where 


4V * iff v fl 2 *i d ? 


H' W/WlVc;^^ 


i|»^ and ^ are the wave functions of the system consisting of the 
ion and the two surface atoms before and after sputtering, respectively. 

is the perturbation (operator) in the Hamiltonian of the 3-body 
system which causes the transition i+f. Equation (20) gives the 
density of final states per unit energy for the normalization volume Q. 


From Eqs. (18) - (20) results the probability for the sputtering 
transition per unit time, 

w 2 (E) = ^ [ft/( 2Trtl) 3 ] 2 |mJ 2) | V(^^~2m) 3/2 (E - E o ) 2 . (21) 


With the assumption w 2 (E) « w^E), one obtains from Eqs. (3), (17) 
and (21) for the relative sputtering probability the approximate 
expression 
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j ; 

i : 


\ ■ 


i 




W S (E) 


a 


iM 


C2) 

if 


8n/2 l^if ; 


(: 


M 


(1) 1 2 m + 2M 


7) 


3/2 


(E - E o )‘ 

J72“ 


( 22 ) 


Since sputtering is a 3-body while reflection is a 2-body interaction, 
the ratio, 

l M if 2> | // Kf^f 2 5 H 21 S/JJ » < 23) 

is essentially the probability w = fi/ft for finding one. surface atom 
within the interaction volume 

fl £ (ft 2 /2mE) 3/2 . , (24) 


the radius of which is of the order of the de Broglie wave length of 
the incident ion. Accordingly, the relative sputtering probability 
in Eq. (22) can be written as 


W S (E) £ 



r (M/m) 2 ,3/2 

l l + 2 (M/m) 


(E - E o )' 


(25) 


The determination of the factor ~ 1> which may be weakly energy 

dependent, requires introduction of an appropriate interaction 
potential and a detailed evaluation of the matrix elements. 
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3. SURFACE SPUTTERING RATE 


In 1) and 2) it has been demonstrated that the statistical method 

[Eq. (13)] and the perturbation approach [Eq. (25)] lead to practically 

the same relative probability for sputtering W (E) . Combining this 

s 

result with Eq. (4) yields for the sputtering rate the formula: 


S(E) a ^ c(E)N 2/3 


r (M/m) 2 .3/2 (E - V 
L 1 + 2(M/ra) J _2 

ft 


(26) 


Since E = 4E s 12 - 35 eV for various metals and E > E for low 
OS ** o 

energy sputtering, Eq. (26) can be simplified to 


S(E) 



a(E Q )N 2/3 


r (M/m) 2 3/2 

L 1 + 2 (M/m) 


< E - V 2 

E o 2 


(27) 


In the considerations under 1), 2), and 3), the effect of the particle 
spin I on the various probabilities has not been included explicitly 
for reasons of a simple notation. In cases where the ion and surface 
atoms have no spin, 1^ =* I g = 0, or the same spin, I 1 = I ^ 0, no 
corrections arise in the expression for W (E) and S(E). In case 
of different spins, I ± ^ I g , W g (E) and S(E) are increased by the 
factor (statistical spin weight) 

S s " 2I g + 1 . (28) 

Equation (27) is exactly of the form of the phenomenological 
expression found by analytically fitting the experimental data.- - -'* 
According to the more general Eq. (25), S(E) reache-o with increasing 
E a plateau-like maximum and decreases then since o(E) decreases at 
sufficiently high energies E. Thus, Eq. (26) appears to be correct 
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not only in the low energy region E > E q , but seems to agree also at 

high energies E » E q qualitatively with the experimental sputtering 

1 - 2 ) 
curves. 

2 2 

It is noted that the particular energy dependence ~(E - E q ) / E 
in Eq. (26) is due to the 3-body interaction which has been assumed to 
be th essential mechanism in low energy surface sputtering. At 
higher ion energies, higher-order many body interactions are energeti- 
cally possible so that not only surface but also an increasing number 
of volume atoms are sputtered (S(E) > 1), 
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4. APPLICATION TO ELECTRICAL DISCHARGE SPUTTERING 


In low pressure discharges, such as glow discharges, practically 

the entire external voltage U Q drops across the cathode sheath X. 

In this region, the discharge ions, which have a thermal velocity 

distribution at the temperature T ■ T fi if the sheath is stable, 

are accelerated to a meat - 1 it r *.ty <v>. The resulting quasi-thermal 

ion beam bombards the ‘ v ..***■. r.V.% ?hlch emits atoms in accordance with 

Eq. (26). Since the ru^hede i'rop is of the order AU ® 20 - 30 volt, 

the unidirectional 'on beam energy is frequently smaller than the 

sputtering threshold, eAU < E q . For this reason, it is essential to 

take into consideration the thermal velocity distribution of the 

ions. Mainly the ions of the tail of the velocity distribution cause 

sputtering of the cathode in this case. 

The cathode sputtering in low pressure discharges is a process 

caused by ions of low energy, E > E q , Accordingly, the number of atoms 

1 2 

sputtered on the average by an ion incident with the energy E - mv 
is given by Eq. (26), 

S(y mv 2 ) = a(|inv 2 - E q ) 2 , (29) 


2/3 H 21 
a = o(E o )N Z/3 -g- 


r (M/m) 2 .3/2 _ -2 

L 1 + 2 (M/m) ^ o 


(30) 


The velocity distribution of the ions, which arrive with the mean 
velocity <v> at the cathode, is assumed to be a Maxwellian in the 
beam system, l.e. 


f (v) = n(m/2irkT) 3 ^ 2 e m ^ V 


<v>) 2 /2kT 


(31) 
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where v designates the Individual velocities of the ions In the 
system of the discharge, and n and T are the density and temperature 
of the ions, respectively. 

The ion beam velocity <v> is determined by the electric field 

+ , 7) 

F in the cathode sheath which is, according to probe measurements,— 

a linear function of the axial discharge coordinate z» 

F(z)=F(l--), F £ 2 U 0 A , (32) 

where X is the extension of the sheath and U q is the external 
voltage. The potential difference across the sheath is 

AU = FX/2 £ U • (33) 

0 


Since collisions involving ions represent a small effect, the magnitude 
of the ion beam velocity at the cathode is in good approximation given 
by 


| <v> | = (2eAU/m) 1/2 « (2eU /m) l/2 


(34) 


The ions within the cone 0 < 0 < tr/2 [0 $ (v, <v>)] which are at 
a distance Az = v cos0 At from the cathode, strike the latter within 
unit time At = 1. Accordingly, the number of atoms sputtered per 
unit surface area and unit time out of the cathode is 


dN _ 

dt " an v 2irkT 


, m .3/2 V* 7 “ m (v - <v>) 2 /2kT 

i 1 J e 


0 0 


*(i mv 2 - E o ) 2 (vcos0)v 2 sin0 d(j> d0 dv 


v = (2E /m) 
o o 


1/2 


(35) 
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by Eqs. (29) - (31). By means of the substitutions, 
cos 0 “ t, dr “ - sin 0 d0 

/iSv “x, dx 13 /fJdv 

and 


(36) 


0 = m/2kT, y => | <v> | , 

x q » /0(2E o /m) 1/2 , y Q = x o - Y , (37) 

Eq. (35) is transformed to 

H - 2 ’ at, <2ST> 3/2B " 2 [ J 0 <E o' <7> - T) E o 2 - J l (E o’ < '" > ' T) ” E ' 1e o 
+ J 2 (E o ,<v>,T)(f ff 1 ) 2 ] (38) 


where 

J n (E o ,<?>,T) - J / x 3+2n e' (x2 + 1,2 - 2Yxt) tdt dx, (39) 


with n = 0,1,2. The T-integration reduces this double integral to 

-v 2 “ 2 

r /r. m \ e t 1+2 n -x , 

J (E q , < v>,T) - — — / x e dx 

x q 

+ ^rJ (y + v) 2(1+n) e _y dy 

1 7 t a. \l+2n -y 2 , 

9 J (y + Y> e J dy 

4y“ y D 


( 40 ) 


4 C 


llpoli a binomial expansion and evaluation of the resulting integrals 
in Eq. (40), there follows; 

-2 , 2 . 2 


*) 


J 0 (E o » < v > ,T) « e" (x o + Y ^ ~ [1 - 'Ky 0 )] 


-2 

8 


- V fCl - 4y 2 ) - 2yy 1 e V ° 


(41) 


“ *2 ty * 2 2k r- 

J 1 (E o ,<v>,T) - 2g- (1 + x 2 ) e~ U ° + Y } ~ Y (5 + 2y‘)[l - #<yj] 


-2 2 
- ^ [(1 - 5y 2 - 8y 4 )-12y 3 y o +(l - 8y 2 )y o -2yy 3 ] e~ y ° 


(42) 


J 2 (E o ,<v>,T) = (2 + 2x 2 + x 4 ) e“^ X o + Y ^ 

+ y (35 + 28y 2 + 4y 4 )tl - <Ky o )l 


+ [(-2 + 14 y 2 + 35y 4 + I2y 6 ) 

+ (35y 2 + 30y 5 )y o + (-2 + I4y 2 + 40y 4 )y 2 
+ 30y 3 y 3 + (-1 + 12y 2 )y 4 + 2y y 3 ] e” y °. (43) 


*) 


In accordance with 


r n -x 4, , -yp v /■, o. v 0 -(l*Hn) n-2m-l 

J x e dx = -e £ (n-l;-2;m)2 'x 

m=0 

+(l~s) (I;2;r)2 ^ +r Vir 4>(x) + const 


r-1 


where (n*,d;m) =* n(n + Id) (n + 2d) . . . [n + (m-l)d] and n ^ 2r-s 
with s = 0 or s » 1. 


4 7 


whore 


4>(y) 


5 jf e" x2 dx » -2- l (-D m y 2m+1 /(2m+l)m! 


/S 


/it m®0 


(44) 


is the error function. Substitution of Eqs. (41) - (43) into Eq. (38) 
leads to the final expression for the number of atoms sputtered per 
unit surface area and unit time from the cathode : 




(45) 


where 


A(E q ,<v>,T) = [1 - 2(1 + x 2 ) x~ 2 + (2 + 2x 2 + x 4 ) x q 4 ] e ( ' X o +Y * 

+ — [1 - (5 + 2y 2 ) x“ 2 + |-(35 + 28y 2 + 4 y 4 ) *“*] [1 - #(y 0 )l 
“3 o 


+ ( (-1 + 4y + 2Yy Q ) 

+ 2 [ ci - 5 y 2 - 8y 4 ) - i2y 3 y Q + d - Qy 2 ) y 2 - 2yy 3 ] x" 2 


+ [ (-2 + 14y 2 + 35y 4 + 12y 6 ) 

+ (35y 3 + 30Y 5 )y Q + (“2 + 14y 2 + 40y 4 )y 2 


+ 30 Y 3 y 3 + (-1 + 12Y 2 )y 4 + 2 Y y 3 ] x^ 4 } e~ y ° 


'46) 


and 


Y = (m<v> 2 /2kT) 1/2 , x 5 (E /kT) 1/2 , 


y Q 5 (E Q /kT) 1/2 - (m<v> 2 /2kT) 1/2 


(47) 
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In Fig, 1, A(y,x 0 ) is shown quantitatively In dependence of the 
2 2 

energy ratios y and x Q . In general, A(y,x q ) increases with in- 

2 2 2 -I 

creasing y and decreasing x q , except in the region y <_ 10 , 

The formula for dN/dt in Eq, (66) holds only at the surface of the 

cathode or within a mean free path from the cathode. In comparing 

this result with experiments, it should be noted that the background 

gas reduces the number of sputtered atoms observed os its pressure 


increases , 
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III. TRANSPORT AND DEPOSITION OF SPUTTERING PRODUCTS 

As long ns the mean free path £ of the sputtered particles is 
very large compared to the distance d between the emitting plane 
and the surrounding system surfaces, the particles travel undeflected 
along straight lines determined by their initial velocity at the point 
of emission. Within this free particle flow, only those system 
surfaces (a) are reached which can be connected with the emitter 
plane through straight lines (Fig. 1), Particles interact, however, 
always more or less weakly independent of how low their concentration 
is since the interaction forces (polarization forces, electric and 
magnetic dipole forces, coulomb forces) have infinite range. For 
this reason, always a few particles will be sufficiently deflected out 
of their initial path so that they can reach system surfaces (b) which 
are not "seen" along a straight line by the emitter (Fig. 2). For 
sputtered particles with a mean free path £ >> d, the deposition on 
surfaces of type (a) can be calculated in first approximation by free 
particle flow, whereas the deposition on surfaces of type (b) has to 
be evaluated by means of a weak interaction diffusion theory. The 
determination of the diffusion coefficient for the sputtered particles 
requires the solution of kinetic equations describing the weak but 
many-partlcle interactions at on the average large distances. The 
analysis of deposition by diffusion on Bystem surfaces of the type 
shown in Fig, 2 leads to a multi-region boundary-value problem with 
mixed boundary conditions. For these reasons, deposition by diffusion 
will not be treated. 
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Z 



FIG. 1: Geometry of emitter (z » 0, 0 £ R <_ R q ) and deposition 

plane (z = d , 0 _< r <_ <*) . 
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Within a cylindrical coordinate system (R, 3, z) , consider an 

emitting surface of radius R ■ R in the plane z = 0. Let the 

o 

system surfaces, for which the deposition is to be determined, be 
located in an arbitrary plane z = d > 0 (Fig. 1) . This general 
deposition problem is solved by calculating the deposition j(r, a) 
in an arbitrary point (r, a, z = d) of the infinite "control plane" 
z = d, 0 < r < ®, 0 <_ a <_ 2 tt. The total deposition on a system surface 
r l £. r 5. r 2 » — a — a 2 * z “ d * s t ^ ,en °b ta:I - ne d by integrating 

j(r,a) over all points (r,a) lying within its boundaries. For ideal 
free particle flow, the depositions on different finite system 
surfaces in the plane z - d do not affect each other or the emitting 
surface. The geometry of Fig. 1 is representative, e.g., for the 
accelerating grid on an ion thruster from which sputtered atoms are 
deposited downstream on system surfaces somewhere in the control plane 
z = d. 

In the plan z = 0, the emitter surface 0 <_ R £ R q , (]<_$<_ 2 tt, 
may emit 

4> = 4>(R, 3) [cm ^ sec (1) 

particles per unit surface and unit time. The rate of deposition at 
the point (r, a) of the control plane z = d due to a differential 
source area do is 

cos 0„ 

dj = *(R,p) £(0^ R,3) 

r l 


do 


( 2 ) 
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where e(0^, R> 3) is the normalized distribution describing the 

emission of the particles in the direction 0^ (Fig. 1) at the point 

(R,8) . It is usually assumed that e(0p R,8) is independent of 

1 ) 

position (R,B) and given by Lambert's cosine law,—' 

e(0^,R,8) = cosO^ir . (3) 

In referring to Fig. 1, it is seen that the following relations hold 
for geometrical reasons, 

9^ «* 0g = 9, cos 0 = d/r^ , 

2 2 2 

s » r + R - 2rR cos cj> , 

_ 2 2 , .2 .. 

r^ =* s + d , <J> * B - a , 


do = RdR dfi 


(4) 


Accordingly, the Impingement rate J(r,ct) in the point (r,a) of the 


control plane z 

J(r,a) - d f f 


d from the entire surface (irR Q ) of the emitter is 


2ir o 


e (9 ,R, B)$ (R,B) RdRdg 


0 [r 2 +R 2 4-d 2 -2rR cos (B-a) J 3 ^ 2 


(5) 


The total deposition per unit time on a system surface bounded by the 
radii r^ and r ^ and the rays and is 


D(r a ) = / / J(r,a) dr da 

1> x ar 

1 r l 


( 6 ) 


Eq. (5) contains as a special case the deposition equation originally 
derived by von Hippel.— ^ Equation (5) can be integrated in closed 
form for cases of practical interest as will be demonstrated next. 
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X) Uniform Emission Source 

In case of a uniform source of sputtered particles which are 
emitted in accordance with Lambert's law, one has [Eqs. (1) and (3)1 s 


$(R,6) - * 0 , 


e(0) = cosO /tt 


(7) 


Hence, 


J(r ,o) 


,2 2ir-« R o 

i • ♦ ; / 

’ ° -c, 0 


R dR difi 


[r 2 +R 2 +d 2 - 2rR cos <fr] 2 


( 8 ) 


by Eq. (5) where $ => 0 - a. The integral is transformed by means 

-1 

of the substitution z = exp(i<|>), d<£ = dz/iz, cos <() = (z + z )/2 as 

di . ,* dz 


2tt-o 
I = / 


2 * $ -1 *> ’ 
-a (p - q cos <j>l |z|=l z ^ p ” ^ z+z 


_ 2 , _2 2 
p = r +R +d , 


2rR 


(9) 


where the point z describes the unit circle as <|) moves from -a 
to 2ir-a. The integrand in the complex z-plane has poles of second 
order at each of the points 

z 1 2 = -[-p + (p 2 -q 2 )^]/q . (10) 


Since p > q > 0, the pole lies outside the unit circle. The 


residue of the pole z^ is 


R 


lim 

z+z. 


4 d r z(z - 2 i>‘ -1 

2 dz , a 2 f ,2 
q L (z-z 1 ) (z-z 2 ) J 


4 *1^2 

2 ( 

q (z r z 2 ) 


( 11 ) 


After replacing z^ and z 2 in accordance with Eq. (10), Cauchy's 
integral theorem yields 
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I ° - i 2triR “ 2Trp/(p 2 -q 2 ) 2 ^ 2 


( 12 ) 


Hence, 


, v ,2 ?° (r 2 +R 2 +d 2 ) dR 2 

J(r,a) * d 0 J 2 ~~2 ~2~2 2 2 3/2 

0 0 [(r Z +R +d ) -4rV] ' 


2 ) 


(13) 


This is an integral of standard type- which is readily evaluated. 
The resulting impingement rate is 


, 2 2 2 
R_ - d - r 

J(r,a) = 1& Q J 1 + ■ " - 2 " ' 2 


2 7 a 2 ,2. 1 fiil 


(14) 


t (R Q + d“ - r )“ + 4r d ]' 


The impingement rate at the point (r,a) of the control plane z = d 
is independent of the azimuthal position a since e and 4> are 
independent of (3 [Eqs. (7) and (9)]. 

2) Parabolic Emission Source. 

For a parabolic source of sputtered particles which are emitted 
in accordance with Lambert’s law, the Eqs. (1) and (3) become 


4(R,8) = 4.(1 - — ,), e(0) = cos8/t r 

V 


. (15) 


Hence 


,2 2ir-a [1-(R 2 /R 2 )1 R dR d<j> 

«*.«) - T ♦. / 


2 2 2 2 
-a [r +R +d - 2rR cos if>] 


(16) 


by Eq. (5) where => 3 - a. Upon substitution of the ^-integral 
evaluated in Eq. (12), Eq. (16) is reduced to 


J(r,a) = h 4 q s 1 + 


R_ 


„ 2 .2 2 
R - d - r 
o 


r 2 . .2 2 2 . 2 , 2^1 

[ (R q + d - r ) + 4r d ] ; 


2 2 2 2 2 
(R + d + r ) R dR 

„ 2 *0 J n r , n 2 2 2.2 2,2 3/2 

R q 0 [ (R + d - r ) + 4r d ] 


4 / 

n y 


( 17 ) 
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2 ) 

The latter Integral is evaluated by standard methods-'. The resulting 
impingement rate at the point (r,a) of the control plane z ■=> d is 


J(r,a) « h$ 0 j 1 + 


„ 4 . . _ 2, ,2 2. . ,,2 , 2.2 

R Q + 2 R q (d “ r ) + (d + r ) 

n 2 , n 2 7.2 272 'f . ijlJi 

R fR + d - r ) + 4r d ] 

o K 0 ' * 
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2 2 

R/ + d 4 

1 + + 2 In f 0 

d 


f 


2 r 2 . ,2 2.2 . 2 . 2 M 

r + [ (R + d - r ) + 4r d ] 


2d 


(18) 


J(r,a) is the same for any azimuth 0 <_ a £ 2m since e and $ are 
independent of 8. The parabolic source distribution [Eq. (15)] is a 
first approximation to the nonuniform emission of atoms from the 
accelerating grid of an ion thruster (presumed tin t: ih density of 
the sputtering charge exchange ions decreases parabolicully with 
increasing R) . 

3) Arbitrary Nonuniform Emission Source. 

In many experimental situations, the source 0 ** HR, 8) is a 
complicated function of both R and 8 with symmetry and boundary 
properties of the form 

*(R 0 ,B) - 0, HR, -8) = -HR, +8) 

d*(0,B)/dR = 0 . (19) 

The emission coefficient is again assumed to be given by Lambert's 
law, 

e(0) = cos0/tt . (20) 

In this case, it is mathematically suitable to expand HR, 3) in a 
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Fourier series, 


CO 00 

♦ ■ ♦„ + n v cos co8 nr* f 

m°0 n=0 o 


(21) 


where 


2tt R o 

I» mn ° / / <I>(R,0)cos mp cos dR dp, 


o 0 0 

o R 
2it o 


* 


o 2irR 


/ / <J> (R, 0) dR d0 


( 22 ) 


o 0 0 


Eq. (21) is a complete expansion for an arbitrary function $»(R,P) 
subject to the conditions in Eq. (19). Insertion of Eqs. (20) - (21) 


into Eq. (15) yields, under consideration of Eq. (14), 

. ( R 2 - d 2 - r 2 

J (r,a) = j ♦„ l + -V 


J 

[(R 2 +d 2 -r 2 ) 2 + 4r 2 d 2 ] 1/2 ) 


+ J(r,a) 
(23) 


where 


0 n 4-l T) 

2 oo oo 2ir-a R o cos m ('( ,+a ) cos ^~ 2 ~ r" RdR d * 

J(r,a) ® ll * f f 2 2 2 S ' 

17 m=0 n=0 mn -a 0 [r +R +d ~2rRcos(J,] Z 

(24) 


By means of substitutions similar to those under a), z = exp(i<)») , . . . , 
cosm(<f> + os) = Re[z m exp(ima)], the <j>-integral in Eq. (24) is transformed 
as 

2tt— ot 


Y = [ cos m((ji+a) dtfr 

m J r „ m2 

-a [p-q cos if»j 
z ( ima 

"4 r ‘ V, f 


m-KL , 
z dz 


z|=l [z^ + 2pq ^ z + 1] 


(25) 


where the integration path is the unit circle in the complex z-plane. 
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The integrand has poles of the second order at each of the points 
z lj2 “ - t-p + (p 2 -q 2 ) 1 ^ 2 ]/q • (26) 


Since p > q > 0, lies outside of the unit circle, the residue 
at the point is 


R 


m+1 . ,2 

Urn d , z lz ~ z l> 


— r - 
1 


z + z l dz l (z-z 1 ) Z (z-z 2 ) 2 

tn r 1 2 , J in 

2 i i* - 3 * / ; 

(Zj-Zg) (Zj-Zp) 


2 , 2 2 , 1/2 

£_ t P “(p -q ) j«> [. 


m 


, 2 2 , 3/2 2 2 J 

(p -q ) p -q 


(27) 


by Eq. (26). According to Eqs. (25) and (27) and the residue theorem 


, 2 2 , 1/2 

P (p q ) m n m 

I » 2r cos mat ] m [■■ ■:■ ■ iH +-~-^].(28) 

q ( p ‘ - q ‘ ) 11 p -q 


m 


Thus, one finds, upon substitution of Eq. (28) into Eq. (24), that 


2 r r . r° ( p(R)“[p 2 (R)"Q 2 (R)1 1/2 


J(r,a) cs 2d \ \ cos ma / 

0 /% IlilS /N 

n=0 . 0 


P(R) 


q(R) 


m 


m 


[p 2 (R)-q 2 (R )) 3/2 p 2 (R)-q 2 (R) 


2n+l_ R n Jn 
cos — — * p" RdR > 

0 


(29) 


where 


p(R) = p(R,r,d) « r 2 +d 2 +R 2 , q(R) = q(R,rd) = 2rR . (30) 


The impingement rate J(r,a) depends both on r and a since 
<5>(R,6) varies with azimuthal position 6. It is seen that the m-th 


61 


Fourier component in Eq, (29) contributes an a- dependence in the 
simple form cos ma. A. two-dimensional emiBBion $(R,B) satisfying 
the conditions in Eq. (19) is, e.g., observed at the accelerating 
grid of ion thrusters if the grid holes are arranged at equal azimuthal 
spacings A<j>(r) along concentric circles r = const < R Q . The 
remaining R-integrals in Eq, (29), in particular those with large m, 
are most conveniently evaluated numerically. 
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IV. KINETIC THEORY OF LOW PRESSURE GAS DISCHARGE 

The theoretical evaluation of electrical discharges by means of 
macroscopic or kinetic equations is one of the classical, unsolved 
problems of plasma physics. To date, no progress has been made in 
calculating the potential distribution In electrical discharges from 
first principles. Although a rigorous theory of the eloctron diode 
with explicit boundary conditions in vacuum based on the Vlasov equation 
is known,— ^ n treatment of the corresponding electron-ion diode in 
vacuum (no ionization) is missing in the literature. In the following, 
a kinetic theory for a low pressure discharge will be formulated. The 
mathematical difficulties resulting from the production of the ions by 
volume ionization (electron-neutral collisions) and the complex func- 
tional boundary conditions for the electron and ion velocity dlstribu- 

i \ 

tions are discussed within the frame of the Vlasov theory,—'' 

In an electrical discharge, the current carriers are generated 

both within the volume of the plasma and at the electrodes. As the 
pressure decreases, the electrons are produced mainly at the cathode 
by secondary processes and thermal emission, whereas the main Bource for 
‘.ons is still volume ionization since emission and production of lono 
at the electrodes is negligible. For this type of low pressure discharge, 
the velocity distributions of the electrons and ions are determined by 
Vlasov equations which are coupled by the Poisson equations for the 
self-consistent electric field and tft sources due to volume ionization. 
The corresponding nonlinear boundary-value problem has functional 
boundary conditions which result from the various electrode processes. 
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In order to reduce the mathematical formalism, a one-dimensional discharge 
geometry with parallel plate electrodes of infinite extension is assumed. 

Consider a low pressure discharge plasma [Debye radius D » 

(V4xr c /kT )*”’* «n] between a plane cathode (x»0) and a plane anode 

SB B 

(x B a) , The mean free paths of the electrons, ions, and neutral atoms 
are assumed to be large compared to the electrode distance, » a 
(so-called "collision-free" plasma). In the one-dlmenBional case, 
the velocity distribution functions of the electrons, f « f(v,x), and 
the ions, F = F(V»x), in the self-consistent field E «= -dO(x)/dx 
are described by the Vlasov equations^ 


9_f , £ d$ 3f_ 
v Dx + m dx 3v 


X a(v,x) 


( 1 ) 


3F e dO DF . 

Dx “ M dx 3V 


S(V.x) 


( 2 ) 


where 


2 +» +* 

=-|* - 4ire / f(v,x) dv - 4ire J F(V,x)dV (3) 

dx -<*> -°° 

is the field-source equation. The electron and ion sources due to 
volume ionization can be reduced to the expressions, 

cr(v,x) « vf(v,x) H[v - , (4) 


E(V,x) » vg (V) / f(v,x) dv , (5) 

—00 


where, 

v o N q < o(v) v > 


( 6 ) 
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g o (V) « (M/2itfeT Q ) 3/2 

-*lMV 2 /feT 

w 

. (7) 


H ( v - V Q ) « 1, 

v > v >0 
- 0 

» 

° 0, 

v < v > 0 
— 0 

* (8) 

is the Heavyside step function. 

In Eqs. (1) - 

(6), e ° -e, and 


e^ = +e, ore the charges, m and M are the masses, and v and V are 
the individual velocities of the electrons and ions, respectively, N q and 
T q are the density and temperature of the neutrals, and a(v) is their 
ionization cross section in dependence of the electron velocity which 
is assumed to be large compared to the heavy particle velocity, 

| v j >> | V | . I is the ionization energy. 

Within the self-consistent field model for long range Coulomb 
interactions, usually (elastic and Inelastic) short-range binary 
interactions are not regarded [X •+• 0 in Eqs. (1) - (2)]. The 
Boltzmann collision integral destroys the "simple" iutegral-functional 
structure of the original Vlasov equation. For this reason, the 
Boltzmann collsion integral for ionization has been approximated in 
Eq, (1) by the discontinuous relaxation expression in Eq. (4). The 
associated ion source in Eq. (2) is then given by Eq. (5) with g Q (V) 
being the Maxwellian in Eq. (7). (The latter presumes that the neutral 
atoms are in thermal equilibrium.) It should be noted that volume 
ionization must be included in the physical model for the discharge 
since otherwise only a trivial solution F = 0 of Eq. (2) would exist. 

On the other hand, the elastic binary interactions of the electrons or 
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ions with the neutrals can be neglected compared to the self-consistent 
long range interactions. 

The boundary conditions for the potential ♦ (*) are given by the 
fixed potentials at the cathode (4> ) and anode ( l t> ) . In contrast to 

O a 

0 and 4> , the boundary conditions for the distributions f(vx) and 
o a 

F(V,x) are not known explicitly. The latter boundary conditions are 
given implicitly by the surface processes at the cathode and anode of 
the low pressure discharge, l.e.s 


i) Reflection of electrons at the cathode (x=*0) and anode (x«a) 

[reflection coefficient: R® «* R e (h mv 2 )] . 

A o 

j) Thermal emission of electrons at the cathode [thermal emission dis- 
tribution: f T = mv 2 )]. 


k) Secondary emission of electrons at the cathode by incident ions 

[probability of emission of an electron of energy h mv 2 by an ion 
of energy MV 2 : S = S(% mv 2 | JjMV 2 )] . 


Z) Neutralization and absorption of the ions at the electrodes [effective 
reflection coefficient: R* “ R* Qi MV 2 )]. 

Accordingly, the boundary conditions for the fields f(v,x), F(V,x), 
and <I>(x) are: 

f(x“0, v > 0) = f T (% mv 2 ) + R^(Js mv 2 ) f(x=0, v < 0) 

+ S(% mv 2 |%MV 2 ) F(x=0, V < 0) , (9) 


f(x=a, v < 0) 
F(x=0, V > 0) 
F(x=a, V < 0) 


= R e (^ mv 2 ) f (x=a, v > 0) 

fi 

« R^(};S MV 2 ) F(x=0, V < 0) £ 0 

» R 1 Os MV 2 ) F (x=a , V > 0) £ 0 
a 


, ( 10 ) 

, ( 11 ) 

, ( 12 ) 
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* 


and 

<t«(x“0) ° $ , iKx=a) = 4> . (13) 

0 u 

Thus, one finds that a low pressure discharge is described by a 

boundary-value problem for nonlinear integro-differential equations 

[Eqs. (1) - (3)] with regular [Eq. (13)] and functional [Eqs. (9) - 

(12)] boundary conditions. For statistical equilibrium and negligible 

2 ) 

field emission, f^, is given by— 

f = B ~0i mv2 + e m " - R e (J$ mv 2 )] H(v) (14) 

1 2irfi 


where T 
metal and 
tunneling , 


is the cathode temperature, e m is the barler energy of the 

£ s 5 ^ is the Fermi energy. Neglecting quantum mechanical 

3) 

the reflection coefficient of the electrodes is— 


R e = , h. mv 2 ^ 1/2 _ ^mv^ 1/2^4 

£ £ 
m m 


(15) 


On the other hand, nearly all the ions are neutralized as they approach 
a metal surface so that their effective reflection coefficient is 

R 1 s 0 . (16) 

2 2 

The secondary electron emission probability S(*s mv \H MV ) has to be 
determined from experimental data since no convincing theory is available. 
The Eqs. (14) - (16) give essentially the correct magnitude of the 
surface effects. A more sophisticated description is not attempted 
since this would render the boundary conditions too complicated. The 
Eqs. (1) - (16) represent probably the most simple theoretical 
formulation of the physical problem under consideration. 
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In spite of the physical simp]/, atlons, the above nonlinear 
boundary-value problem turns out to be extremely difficult. This is 
mainly due to the source terms cr(v,x) and E(V,x) in Eqs . (1) - (2). 
If one treats volume ionization as a small perturbation, then Eqs. (1) 
(2) reduce to the ideal Vlasov equations in first approximation (X 0) 


9f ,e d4> 9f n 

v — *r — — — — a n 

9x m dx 9v 


9F e dO 9F 
V 9x “ M dx 9V " 0 


(17) 


which have generalized functional similarity solutions of the form 
f = P[|rav 2 - e<I>(x)], F = Q[| MV 2 + e<Kx)] . (18) 

P and Q are arbitrary discontinuous functions [because of the 
hyperbolic nature of Eq. (17)] of the functional argument indicated 
which have to be determined in such a way that the boundary conditions 
in Eqs. (9) - (12) are satisfied. This perturbation approach does, 

•t 

however, not work in absence of volume ionization (X=0) since 
F(v,x) = Q = 0 because of the homogeneity of the boundary conditions 
for F(v,x) [Eq. (12)]. For this reason, a more general functional 
similarity solution has to be derived for the complete Eqs. (1) - (2) 
with volume ionization (X=l) . This could not be accomplished, however, 
in spite of a considerable effort in time. It is hoped that the 
solution of this gas discharge problem can be reported at a later date. 
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